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Abstract 

A general framework is presented to discuss the approximate solutions of an evo- 
lution equation in a Banach space, with a linear part generating a semigroup and a 
sufficiently smooth nonlinear part. A theorem is presented, allowing to infer from an 
approximate solution the existence of an exact solution. According to this theorem, 
the interval of existence of the exact solution and the distance of the latter from the 
approximate solution can be evaluated solving a one-dimensional "control" integral 
equation, where the unknown gives a bound on the previous distance as a function of 
time. For example, the control equation can be applied to the approximation meth- 
ods based on the reduction of the evolution equation to finite-dimensional manifolds: 
among them, the Galerkin method is discussed in detail. To illustrate this framework, 
the nonlinear heat equation is considered. In this case the control equation is used to 
evaluate the error of the Galerkin approximation; depending on the initial datum, this 
approach either grants global existence of the solution or gives fairly accurate bounds 
on the blow up time. 
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1 Introduction. 



In this paper we consider, within a Banach space F, a Volterra integral equation 

<p(t)=U(t-t )f + I dsU(t-s)P(<p(s),s) , (1.1) 

J t 

for an unknown function ip from a real interval to F. Here /o S F, IA is a linear semigroup 
on F and is a locally Lipschitz nonlinear map from an open set of F x R to F. If IA is the 
semigroup generated by a linear operator A : DomA CF-»F, under minimal technical 
conditions the above Volterra equation is equivalent to a Cauchy problem 

<p(t) =Atp(t)+V{<p(t),t), <p(to)=f , C:=d/dt). (1.2) 

To standardize the language, problems (jl.lj) . (|1.2j) are defined precisely in SectEl local 
existence and uniqueness of their solutions are well known. 

The aim of this paper is to discuss the approximate solutions of (|1,1[) . In the most general 
sense, an approximate solution is simply a continuous map t \— * <p ap (t) which can be 
inserted in the r.h.s. of (|1.1|) . i.e., such that graph p ap C DornP. For any such map, we 
can define the integral error as the difference between the two sides of If ip ap is a bit 

more regular, the integral error is determined by the differential and datum errors which 
are, respectively, the differences between the two sides in the differential equation and in 
the initial condition of (|1.2j) . 

All the above concepts are formalized in SectEl Here, we also present a general statement 
(Prop 1231 which can be applied to an approximate solution t t— > <p ap {t) to infer the ex- 
istence of an exact solution ip on an appropriate time interval, and also to estimate the 
difference cp(t) — (p ap (t). The essential character in Prop l3.4l is an integral control inequal- 
ity, depending on the available estimators for the integral error of ip ap and for the growth 
of V away from the graph of ip ap . 

The unknown in the control inequality is a real, nonnegative function t \— ► R(t); if a solution 
R is found to exist on a time interval [to, t\ | (i.e., either [to, ti] or [to, t\)), then it is granted 
that (jl.lj) possesses an exact solution ip : [to,t±\ — > F, and that \\<p(t) — p ap (t)\\ < R(t). 
In typical cases, a solution of the previous integral inequality can be constructed solving 
an ordinary differential equation for R, that we call as well the control equation. In this 
way, the problem of giving estimates on the existence time for (|1.1|) and on its exact 
solution p, living in F which is typically of infinite dimension, is reduced to the analysis 
of a one-dimensional ODE. 

Prop l3.4l can be regarded as a general formulation of many statements about specific 
evolutionary problems, often encountered in the literature. From this viewpoint, the 
content of this Proposition is not at all surprising: however, the technique we use to prove 
it is essentially different from the arguments often employed in related situations. The 
standard way of thinking would suggest to prove Prop 13.41 in two steps: a) derive (via some 
nonlinear Gronwall lemma [T6j) an a priori bound \\p(t) — cp ap (t)\\ < R(t), holding until 
p(t) exists; b) show that nonexistence of p on the whole interval [to 5 * 1 1 would contradict 
the previous bound: this argument is called the "continuation principle" in |19| . 
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On the contrary, the proof we propose (in SectEJ is very direct, and shows that 92 can 
be constructed on the whole [*q , tx I by a convergent Peano-Picard iteration, applying 
repeatedly the Volterra integral operator to the approximate solution ip av . The control 
inequality ensures the invariance under the Volterra operator of the space of functions 
with distance < R(t) from <p ap (t) on [£o,ii|; the confinement to this domain of all iterates 
of (p ap , and the local Lispchitz nature of V, allow to prove their convergence to a function 
ip, also distant less than R from ip ap . 

As a first, very simple illustration of PropEHJ in Sect we apply the control equation to 
the approximate solution tp ap {t) := 0. In spite of the trivial choice for (p ap , the control 
equation gives useful information on the interval of existence and on the growth of the 
exact solution (p, depending on the norm ||/o|| of the initial datum. The accuracy of these 
predictions is tested on an example, concerning the (one-dimensional) wave equation with 
polynomial nonlinearity. 

A second, more refined application of the control equation is proposed in SectE] for the 
Galerkin scheme (and similar approaches). In the conventional formulation, the Galerkin 
method is an algorithm to construct approximate solutions 1 1— > <p ap {t) of with values 
in a finite-dimensional submanifold of F. In this Section, the standard evolution equations 
for the coordinates of tp ap (t) in the Galerkin submanifold are coupled with the control 
equation for R(t); in this way, a finite-dimensional system of ODE's gives simultaneously 
the Galerkin approximate solution <p ap , an interval [io>^i| on which the exact solution <p> 
of is granted to exist and an upper bound for ||</? ap (i) — ip(t)\\ on this interval. 
In Sect [71 all the previous results are applied to a nonlinear heat equation, working for 
simplicity in one space dimension (with a spatial coordinate x G (0, tt)). In this case, the 
Cauchy problem (|1.2|) (with initial time to := 0) is, symbolically, 

<p{x,t) = <p xx (x,t) + ip(x,t) p , cp(x,Q) = f (x) (1.3) 

with p G {2,3,4, ...}, to be discussed in the Sobolev space F := Hq(0,tt). The implemen- 
tation of the general framework in the present case with polynomial nonlinearity requires 
accurate information on the pointwise product of functions in Hq (0, 7r); in particular, pre- 
cise estimates are needed for the norm \\fh\\ when /, h are in this space (see the AppendixlAl 
about this, and |llj for more general information about multiplication in Sobolev spaces). 
To exemplify some general facts about (|1.3|) . in the same Section we consider the initial 
datum fo(x) := sinx. If the (nonnegative) constant A is below a critical value, the 

control equation for the zero approximate solution suffices to prove existence of a globally 
defined solution ip : [0, +00) — > F of (|1.3|) . For larger A, the same control equation 
gives a finite lower bound for the existence time of the solution (p. These conclusions are 
complementary to the ones arising from a known "blow up" theorem of Kaplan for the 
nonlinear heat equation (see [Sj; a review is given in the Appendix [Bj) . When Kaplan' s 
theorem is applied to ()1.3|) with the previous datum, for sufficiently large A it predicts a 
finite, explicitly determined upper bound on the existence time of the solution. 
Again in SectJTJ we add to the above facts the information arising from application of the 
control equation to the Galerkin scheme; the chosen Galerkin submanifold is the linear 
span of finitely many elements in the Fourier basis. As an example, we consider the 
Galerkin differential equations for two modes, coupled with the control equation for R, 
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with p = 2 and the previous f$. This system in three unknown real functions can be 
easily treated by any package for the numerical solution of ODE's; the results obtained 
by the MATHEMATICA package, for several values of A, are presented with some detail. 
Among other things, the Galerkin approach with the control equation allows to increase 
the critical value of A below which global existence is granted for 1)1.3)) : for A above the 
new critical value, a better lower bound for the existence time is derived. If A is fairly 
large, the new lower bound is close to the Kaplan upper bound, which yields an uncertainty 
between 20% and 30% on the existence time of the exact solution. Also, the upper bound 
R(t) on \\<p(t) — <p ap (t)\\ is fairly small in comparison with ||y op (t)|| for non large t. 
To some extent, it is surprising that a fairly good accuracy can be obtained combining the 
control equation with a Galerkin scheme in two modes only. These outcomes encourage us 
to hope that the same method would give nontrivial information on the Cauchy problem 
for the equations of fluid dynamics, whose Galerkin approximations in few modes give 
rise, among others, to the widely studied Lorentz model |Hj [IB] . 



2 Preliminaries. 



Throughout the paper, F denotes a real or complex Banach space with norm || || and 
elements /, /o, /i, h, ... . We write B(/o, p) for the open ball in F of center /q and radius 
p (if p = +oo, this means the whole F). Let us be given a linear operator 



A : DomA C F 



(2.1) 



with domain a linear subspace of F; whenever we speak of a continuous map from/to 
DomA, we always refer to the topology of the graph norm ||/||.4 := ||/|| + \\Af\\ (as well 
known, DomA is complete in this norm if and only if A is closed). We denote with C(F) 
the Banach space of bounded linear operators of (the whole) F into itself. 
We always write [io,ti| for a real interval of the form [to,ti] or [to,ti) (always intending 
to < t±) in the second case, ti can be +oo). If if) : [*o? *i I ~> F, the graph of this function 
and the tube around tp of any radius p G (0, +oo] are 



graph ip := { (ip(t),t) | t G [i ,*l| } C F x R . 

Tty,p) :={(/,*) eFx [t 0) h\ | \\f-m\\<P} 

(the latter is the whole F x [to,ti|, if p = +oo; it becomes B(/o,p) x [to,ti| 
const. = f ). 

Linear semigroups on F. This name indicates maps U such that 

U : [0,+oo) -»■ £(F) , U(t + s)=U{t)U{s), U(0) = 1 F . 
The generator of a linear semigroup IA is the linear operator 

A : DomA CF^F, / h+ Af , 

d 



(2.2) 

(2.3) 
if V(t) = 



(2.4) 



DomA := {/ G F 



d 
It 



t=o 



U(t)f exists} , Af :-- 



t=o 



U{t)f 



(2.5) 
(2.6) 
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(with (d/dt)t=o denoting the right derivative). 

A linear semigroup U in F is strongly continuous if for all / G F the map [0, +00) — > F, 
t 1 ^ U{t)f is continuous. In this case (see, e.g., [2] ^2]), the map (/, t) 1— > U{t)f is jointly 
continuous, the generator .4 is densely defined in F and closed, and a) b) hold: 

a) A determines U. For all /o G DomA, the function 1 1— ► U(t)fo is the unique function 93 
such that 

93 G C([0, +00), DomA) n C 1 ([0, +00), F), 0(t) = ,4</>(i) for all t, cp(0) = f ; (2.7) 

b) for any function ?/> G C([io, DomA) D C 1 ([t , £i|, F) and t in this interval, it is 

^(t) = W(i - t )ip(t ) + / dsU{t-s) tp(s) - Atp(s) 

Jt L 

(here and in the sequel, the dot indicates the derivative). A linear semigroup IA is uniformly 
continuous if the map IA is continuous from [0, +00) to /2(F) with the standard operator 
norm (uniformly continuous semigroup); this happens if and only if IA has generator A G 
C(iF), and gives a trivial example of strongly continuous semigroup (extendable to t < 0). 



2.1 Definition. An estimator for a linear semigroup U on F is a continuous function 
u : [0, +00) — > [0, +00) such that, for all f G F and t G [0, +00), 

[|W(t)/[|<«(f)||/||. o (2.9) 

Each strongly continuous linear semigroup admits an estimator of the form u(t) = Ue~ Bt , 
where U > 1 and B are real constants: see [H]. 

Lipschitz maps. If C, T> are subsets of a topological vector space, we say that C is a 
strict subset of T>, and write C <s T>, if C is bounded and C <Z T>, the symbol — denoting 
the closure. Now, let us be given a (possibly nonlinear) map, with open domain, 

V : DomV cFxR->F, (f,t)^V{f,t). (2.10) 



2.2 Definition. We say that V is Lipschitz at fixed time (or, respectively, Lipschitz) on 
the strict subsets of its domain if, for every C <s DomV, there is a nonnegative constant 
L = L(C) (or, resp., a pair of nonnegative constants L = L(C), M = M(C) ) such that 

\\V{f,t)-V{f,t)\\ < L\\f - f'W for (f,t), (f',t) G C; (2.11) 
\\V(f, t) - V(f, t')\\ <L\\f-f'\\+M\t- t'\ for (/, t), (/', t') G C . o (2.12) 

Of course (|2,12|) implies (|2.11|) and the continuity of V . 

An example. Some applications presented in the sequel rely on a map V of the form 
DomV = F x A (A C R an open interval), V(f, t) := T(f, /, t) , (2.13) 
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where 

y : x p F xA-tF (p E {1,2, ...}) , (f u ..., f p , t) -+ 3>(/i, f p , t) (2.14) 
is R-linear in each argument fx, f p ; it is also assumed that 

mil, ~,fp,t)\\<p® \\fi\\-\\fp\\ , (2.i5) 

P(f 1 ,...,f p ,t)-T(f 1 ,...,f p ,t')\\<Q(t,t') \\h\\...\\f P \\ 

for all fi,...,f p £ F and t,t' G A, where P : A -» [0, +oo) and Q : A x A -> [0, +oo) 
are continuous functions. It is finally required that, for each fig A, there is a constant 
M = M{B) such that 

Q(t,t') <M\t-t'\ fovt,t'eB. o (2.16) 



2.3 Proposition. With the previous assumptions, for all f,f'£F and t, t' E A it is 

\\v(f,t)-v(f,t')\\ <P(t)J2( P j ) \\f'\\ p ~ j \\f-f'\\ j + Q(t,t')\\fT • (2-17) 

Proof. Setting for convenience f\ := /', f2'- = f~ f we can write 

r(f, t) - vw, t) = nh + h, h + f2,t) - nfu fut)= (2.i8) 

p 

= E E A pi :={(h,. ..,l p ) e {1,2}P \ l s = 2 for j values of s } . 

j=i {h,...,ip)eA P j 

From the first inequality (|2.15|) . we infer 

\wu,t)-v{f,t)\\<p{t)j^( p . \ \\r\r j \\f-f\\ j , (2-i9) 

j=i 

because A p j has cardinality ^ ^ j . Finally, the second assumption (|2,15|) gives 

\W(f,t)-V{f,t')\\ <Q(M')II/T (2-20) 

and (|2.19|) ()2.2(J|) . with the triangular inequality, yield the thesis (|2.17|) . o 
Ea. (|2.17|) will be frequently used in the sequel; together with (|2.16|) . it implies 

2.4 Corollary. The map V is Lipschitz on the strict subsets of Fx A. o 

General formulation of the Volterra and Cauchy problems. We define formally 
both problems, and review their relations. 
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2.5 Definition. Let us be given: 

i) a strongly continuous linear semigroup U on F; 

ii) a continuous map V : DomV C F x R —* F, with open domain; 
Hi) a pair (/o,to) £ DomV. 

The Volterra problem related to U and V with datum /o at time to ^ s the following one: 
Find ip £ C([to,ti|, F) such that graph ip C DomV and 

<p(t) =U(t -t )/o+ / dsU(t-s)V(ip(s),s) for all t e [t ,h\ . (2.21) 
■/to 

2.6 Definition. Consider: 

i) a linear operator A : DomA C F-> F; 

a continuous map V as in the previous definition; 
Hi) a pair (/o,to) £ DomV such that /o 6 DomA. 

The Cauchy problem corresponding to A, V with datum /o at time to is t/ie following one: 
Find ip G C([to, -Dom^4) n C 1 ([to, ti|, -F) such that graph cp C DomV and 

<p(t)=A<p(t) + V(<p(t),t) foraJlte [to, ti| , <p(t ) = fo ■ (2.22) 

2.7 Proposition. Let A, V,fo,to be as in De f\2.fk and further assume A to be the 
generator of a strongly continuous linear semigroup U. Then: 

i) a solution cp of the Cauchy problem 12. 2$\) is also solution of the Volterra problem 12. 21]) : 

ii) as a partial converse, a solution ip of the Volterra problem 12.21]) is a solution of the 
Cauchy problem 12. 22]) in either of these situations: a) F is reflexive and V is Lispchitz 
on the strict subsets of its domain; (5) (trivial case) A € £{F), no further assumptions on 
F and V. 

Proof. It is essentially based on (|2.7I2.8JI : see j2j (conditions a) (5) in item ii) ensure a 
solution if £ C([t ,ti|,F) of (l2~2H to be in C 1 ([to, ti|, F) n C([t , h\, DomA)). o 

In particular, the operator ^4 := is the generator of the identity semigroup U(t) = lp for 
all t. With this remark, the framework of this paper applies to any ODE cp(t) = V(^p(t),t) 
in a Banach space, also including the finite dimensional cases F = R m or C m . 

2.8 Proposition. Consider the Volterra problem \2.21\) . where li is a strongly continu- 
ous linear semigroup, and V is continuous and Lipschitz at fixed time on the strict subsets 
of its domain; then i) ii) hold. 

i) Problem 12.21]) has a solution. 

ii) If ip : [to, ti| - * F and ip' : [to^il ~~ > ^ are two solutions, it is 

<p(t) = <f/(t) forte [t , ti| n [t , t[\ . (2.23) 
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Proof, ii) We consider any £2 in the intersection of the domains. Subtracting Eq. (|2.2H) 
for ip from the analogous equation for ip' , and taking the norm, we obtain 



\\<p(t)-</(t)\\< dsu(t-s)\\V^(s),s)-V^'(s),s)\\<UL \\<p(s) - ip'{s)\\ (2.24) 
Jto Jto 

for each t £ [to, £2]- Here: it is any estimator for U; U := max se [ fo ^ 2 ] u(s); L is a constant 
fulfilling the Lipschitz condition (|2.11f) for P on the set C := graph (99' \ [to, *2])Ugraph {ip f 
[£0,^2]) (this C is a strict subset of DorriP). Eq. ()2.24jl and the classical Gronwall Lemma 
[nil imply ||v?(t) - ¥>'(i)|| = for all t E [i ,t 2 ]. 

i) Eq. (|2.21|) is the fixed point problem for a Volterra type integral operator, and a solution 
can be constructed by standard Peano-Picard iteration, starting from the function tpo(t) := 
const. := /o; see, e.g., j2j. o 

From our viewpoint, the previously mentioned argument for local existence is a particular 
case of a more general statement, allowing to construct a solution of Q2.21JI by a Peano- 
Picard iteration with starting point any approximate solution (of sufficiently small error); 
all this will be discussed in the next Section. 

Of course, Prop |2.7l allows to transfer the statements on uniqueness and existence from 
(|2.21j) to (|2.22jl . Let us call maximal a solution ip of (|2.21|) or (|2.22|) which has no proper 
extension. If one can grant the existence of a solution on a sufficiently small interval, 
and the coincidence of two solutions on the intersection of their domains, it follows that 
a unique maximal solution <p exists and any other solution is a proper restriction of the 
maximal one. Furthermore, if local existence is granted for arbitray data, the domain of 
the maximal solution ip with a given datum has the form [to,i9) (otherwise, <p> could be 
extended taking its value at 1? as a new initial datum). 



3 Approximate solutions. Statements of the main results. 

We consider a strongly continuous linear semigroup U on the Banach space F, and a 
continuous function V : DorriP C F x R — ► F with open domain. We are interested in the 
Volterra problem (|2.21j) . for a given pair (fo,to) £ DorriP. 

3.1 Definition. By an approximate solution of problem \2.21\) . we mean any continuous 
function (p ap : [to > ill — * F, such that graph (p ap C DorriP. Given any such function, we 
stipulate the following: 
i ) the integral error of <p ap is the function 

E(p ap ):te [t , ti| ^E(cp ap )(t) :=p ap (t)-U(t-t )fo- [ dsU(t-s)P(<p ap (s),s) ; (3.1) 

an integral error estimator for ip ap is a continuous function £ : [to,ti| — > [0, +00) such 
that, for all t in this interval, 

\\E(cp ap )(t)\\<£(t) . (3.2) 
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ii) The datum error for (p ap is the difference 

d(<p ap ) ■= <P aP (to) - fo ; (3.3) 

a datum error estimator for ip ap is a nonnegative real number 5 such that 

\\d( Vap )\\<6. (3.4) 

in) If A is the generator ofU and ip ap £ C([to, t±\, DomA) n C 1 ([to, h\, F), the differential 
error of (p ap is the function 

e(<p ap ) : t £ [t , h\ i-> e(ip ap )(t) := ip ap (t) - A(p ap (t) - V(cp ap (t),t) ; (3.5) 

a differential error estimator for (p ap is a continuous function e : [io^il ~~ > [0, +oo) such 
that, for t in this interval, 

\\e(<p ap )(t)\\ <e(t) . o (3.6) 

Of course, ip ap is a solution of the Volterra (resp., Cauchy) problem iff E(ip ap ) = (resp., 
d((p ap ) = and e(ip ap ) = 0). 

3.2 Lemma. Let (p ap : [to^il F be an approximate solution of the Volterra problem 
\2.21\) . and assume the regularity conditions in item Hi) of the previous Definition. Then, 
the integral error of ip ap is related to the datum and differential errors by 

E{ip ap ){t) =U(t-t )d(<p ap )+ I dsU(t-s)e{<p ap )(s) . (3.7) 

J t Q 

IfU, d(ip ap ), e{ip ap ) have estimators u,5,e, then E((p ap ) has the estimator 

£(t) := u(t - t ) 5 + / ds u(t - s)e(s) for all t e [t ,h\. (3.8) 

Proof. Ea. (|3.7|) follows applying the definitions of E(ip ap ), d(ip ap ), e(ip ap ) and the identity 
with V := 9V- Given (j?TT|) . the bound (|3~£|) on \\E(tp ap )(t)\\ is evident. o 

Remark. The estimator £ defined by (|3,8|) is useful, because in many cases it can be easily 
computed. However, in peculiar situations involving oscillating functions, this estimator 
can be rough. For example, consider the semigroup U{t) := lp for all t, with generator 
^4 = and estimator u(t) := 1. Let us choose ip ap (t) ■= fo for all t, so that d(tp ap ) = 0, 
e(ip ap )(t) = -V(f ,t) and E{ Va3 ){t) = - / t * ds V(f ,s). Suppose V{f ,t) = g e wt , with 
lo £ (0, +oo) and go a vector of the (complex) space F; then, the best estimators for 
d{(p ap ) and e(ip ap ) are, respectively, 5 = and e(t) = \\go\\. Correspondingly, Eq. (|3.8j) gives 
the integral error estimator £{t) = \\go\\{~t — to); on the other hand, it is found by direct 
computation that E(ip ap )(t) = igo[e lwt — e lujt °)/u>; thus ||£J(^ ap )(f)|| is a bounded function 
of t, whereas the estimator £ grows linearly. 

Similar drawbacks of the estimator (|3.8|) in the presence of oscillatory functions are met 
(even for F = R m ) if one considers a differential equation with fast periodic variables and 
the approximate solutions which arise from averaging methods jSJ. o 
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To formulate the main theorem on approximate solutions, we need one more notion de- 
scribing the growth of V away from a function ip : [to, t\\ — ► F, such that graph ip C DomV. 

3.3 Definition. A growth estimator for V from ip (if it exists) is a continuous function 

£ : [0,p) x [t , ti| [0,+oo) , (r,t)^£(r,t) (3.9) 

such that: 

i) p£ (0,+oo] and T{ip,p) C DomV (see Eq. fgjfy ); 

ii) £ is nondecreasing in the first variable: £(r,t) < £(r',t) for r < r'; 
in) for all (f,t) £ T(ip,p), it is 

\\v(f,t)-rm),t)\\ <m-m\\,t) . <> (3.10) 

Remarks, a) From a continuous function £' : [0, p) x [to, ti| — ► [0, +00) fulfilling i) iii) but 
not ii), we can construct the function £(r,t) := max r / g [ ,r] ^'( r ';t), which also fulfils ii). 
b) A function V which is Lipschitz at fixed time on a tube around ip possesses on it a 
growth estimator linear in r. Less trivial estimators appear if DomV = F x A, A a real 
interval, and one wishes to estimate the growth of V from ip on the whole product space 
F x [to ; ^i| (= on a tube of infinite radius). For instance, consider a map V as in the 
Example of page EJ the growth of V from any : [to, t\ \ C A — > F admits the estimator 

£{r,t):=P{t)j^( P .\u(t)\r 3 rl (3.11) 
j=i 

(r £ [0,+oo), t £ [t , t^). To find this, apply En. dTTTll with /' = ip(t) and t' = t. o 
We come to the main theorem of this Section: the proof will be given in Sect|l] 

3.4 Proposition. Let us be given a Volterra problem where: U is a strongly 
continuous linear semigroup; V : DomV C F x R — > R is continuous and Lipschitz at 
fixed time on the strict subsets of its open domain (Def. \2.2\) . Assume that: 

i) u is an estimator for U; 

ii) ip ap : [to,t\\ — > F is an approximate solution of h2.21\) . and £ : [to,t\\ — > [0, +00) is an 
estimator for the integral error E((p ap ); 

iii) £ : [0, p) x [to>ti| - * [0, +00) is a growth estimator for V from ip ap (p G (0, +oo]j. 
Consider the following problem: 

Find R G C([t ,*l|, [0, p)) such that 
£(t)+ f ds u(t- s) £{R(s),s) < R{t) for t £ [to,h\ . (3.12) 

Jto 

Lf \3.12\) has a solution R with domain [io,ti|, then \2.21\) has a solution ip with the same 
domain, and for all t therein it is 

\\<p(t) - <p ap (t)\\ <R(t) ■ (3.13) 
The solution ip can be constructed by a Peano-Picard iteration, starting from ip ap . o 
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3.5 Definition. Eg. hH.lfy) will be referred to as the control inequality. o 

Remarks, i) The function R is required to exist on the same domain [io,^i| of (p ap . In 
many applications, one starts with an approximate solution on a domain [to? ^2 1 an d then 
finds 1)3.12(1 to a have a solution R on a domain [to, t\\ C [to ? *2 1? °f course, in this case the 
previous Proposition must be applied to ip ap \ [to,ti\. 

ii) As anticipated, the argument we will employ to prove Prop l3.4l is different from the "con- 
tinuation principle" mentioned in the Introduction. Instead of using a Gronwall Lemma 
plus a reductio ad absurdum, we will prove the existence of tp on [io>il|> and the bound 
(|3.12j) . by a constructive Peano-Picard iteration; the convergence of this iteration on the 
whole [£o)ti| has some theoretical interest by itself. Furthermore, this approach overcomes 
some technicalities required by the application of nonlinear Gronwall lemmas (the analysis 
of the associated integral equation, and the necessity to determine the greatest solution 
when uniqueness fails jl6j). 

Apart from the general concept of approximate solution employed here, the idea to prove 
existence for an ODE / = V(f, t) by the Peano-Picard method, under conditions of non- 
linear growth for V of more global type than the Lipschitz property can be ascribed to 
Caratheodory pQ, and was developed in ^7] [T^] . 

iii) Of course, we can accept as a solution of ()3.12j) an R fulfilling the equation 

£(t)+ I ds u(t- s)£(R(s),s) = R(t) for t e [t ,h\ , (3.14) 
Jto 

hereafter referred to as the control integral equation. (The existence of such an R on a suf- 
ficiently short interval is granted by standard compactness arguments, see [T|I|. Uniqueness 
can be proved under supplementary assumptions of Lipschitz kind for I). o 

Let us exploit a typical case, where the control integral equation (|3.14j) is equivalent to a 
Cauchy problem. To this purpose, assume that 

u {t) = Ue~ Bt (U>l,Be R), £(t) = Ue- B ^ ta) 5 + U I ds e~ B{t ~ s) e(s) (3.15) 

Jt 

for some constant 5 > and some continuous function e : [to,*i| - ► [0, +oo) (for example, 
the estimator £ derived from Eq . 1)3.8)1 has the above form). Then, multiplying by e B ^~ to ^ 
we see that Eg . 1)3.14)1 is equivalent to 

US + U I ds e B ( s - to h{s) +U I ds e B ^ s - to k(R{s), s) = e B{t -^R{t) . (3.16) 

Jto Jto 

Any solution R of (j3.16)) is clearly C . By derivation in t of this equation, and evaluation 
of the same at t = to , we get 

3.6 Proposition. If u and £ are as in \3.15\) . Eq. \3.lJ\) is equivalent to the problem 

R(t) = Ue(t) + U£(R(t),t) - BR(t) , R(t ) = U5 , (3.17) 

for an unknown function R £ C 1 ([to, h\, [0, p)) (the terms control problem, or control 
equation will be employed as well, for Jff.i7| ) or the differential equation therein). o 
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4 Proof of PropJSm 

We present in detail the argument in the case of a compact interval [io,ti]- In this case, 
we use the space C([to,ti],F), regarded as a Banach space with the usual sup norm 
H^ll := max tg [ to tl ] The case when cp ap , R, etc. are defined on [to,ti) (with t\ 

possibly infinite) is treated in a similar way, using C([to,ti),F) with the topology of 
uniform convergence on all compact subintervals [to,r] C [to,t{) Sticking from now 
on to the case [to>*i]) we introduce the objects 

q : = max R(t) ; T(<p ap , g) = {(/, t) £ F x [t , h] \ \\f - <p ap (t)\\ < 0\ i t 4 - 1 ) 
te[t ,ti] 

D := {tp e C([t ,h],F) | \\i/)(t) — Pa P (t)\\ < R(t) for t € [to, t\] } ; (4.2) 

then g < p, and T(ip ap , g) C T((p ap , p) is a strict subset of DomV. D is a closed subset of 
C([t ,ti],F) (containing Lp ap ) and )/>£3=> graph ip C T(ip ap , g). 

4.1 Definition. We put 

J:Q^C{[to,h],*) > (4-3) 

JW(t) :=U(t-t )f + [ dsU(t-s)V(ij(s),s) Vt€[to,ti]. o 
Of course, we have 

4.2 Lemma, ip S 2) solves the Volterra problem \2.21]) if and only if (p = J{^p) ■ 

4.3 Lemma. There is a constant A > suc/i i/iai, /or a// G 23 > 

W)(*)ll < A /ort€ [t ,ti] • (4.4) 

Thus — J{ip')\\ < A(ti — to)!!^ ~~ ^'IIj which implies the continuity of J . 

Proof. V is Lipschitz at fixed time on the strict subsets of its domain, so there is a 
constant L > such that \\r(f,t)-V(f',t)\\ < L\\f - f'\\ for (f,t), (f',t) e T(p ap ,g). If 
u is an estimator for the semigroup U and £7 := maxjg^ ^] it(i), we see that Ea. (|4.4|) is 
fulfilled with A := UL; the remaining statements are trivial. o 



4.4 Lemma. J(D) C £>. 

1 This complete, locally convex topology on C([to,£i),F) is denned by the seminorms (|| || T ) rg [ to tl j 
where |^|| T := max 16 [ l0)T ] ||i/>(t)||. To adapt the proof to this case, the objects g, T(<p a p,g), A, etc. 
appearing in the sequel must be replaced by families of objects g T , T(<^ ap \ [to,r], g T ), A T , etc., one for 
each r; the definition of D is simply rephrased using \to,t{). The Peano-Picard iteration converges in all 
the seminorms || || T . 
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Proof. Let ip S D. For all t G [io^i] the definitions of J and of the error E(ip ap ), with 
the properties of £,u,£, imply 

J(i/i)(t) - <p ap (t) = -E{ Vap )(t) + [ ds U(t - s) [ Vty{s), s) - V(<p ap (s), s) } , (4.5) 

J t 

\\J(m) - <Pap(t)\\ < £{t) + / ds u(t - s) i{U(s) - <pM\\,a) . (4.6) 

Jt 

On the other hand, || — </j ap (s)|| < R(s) which imphes £(\\ i/>(s) — (p ap (s)\\, s) < £(R(s), s); 
inserting this into (|4.6|) . and using the control inequality (j3.12|) for R, we conclude 

\\J(i;)(t)-<p ap (t)\\ < R(t) , i-e-, JU>) G 3). o (4.7) 

The invariance of D under J is a central result; with the previously shown properties of 
i7, it allows to set up the Peano-Picard iteration and get ultimately a fixed point. 

4.5 Definition. (<pk) {k £ N) is the sequence of functions in D, defined recursively by 

¥>o ■= ¥a P , ■= J((fk-i) (k > 1) . o (4.8) 

4.6 Lemma. For all k £ N and t £ [to,ti], it is 

y k+1 (t)-Mt)\\<^ Ak(t ~ to)k (4.9) 

where A is the constant of Eg. {4-4\) an d S := max tg r t0)tl i S (t). So, 

^Jc ^ \k 

[|yfc+i-y>fcll<S fc , ■ (4.10) 

Proof. Ea. (|4.1U|) is an obvious consequence of (|4.9|) . We will prove (|4.9|) by recursion, 
indicating with a subscript ^ the thesis at a specified order. 

We have cpx(t) - <p (t) = J(ip ap )(t) - tp ap {t) = -E(ip ap )(t) by the definition of E(ip ap ), 
whence ||yi(t) — ¥>o(£)|| < £(t) < S; this gives (|4.9|) n. For each k > 0, we have 

||Vfc+2(*) " ¥>fc+l(*)ll = \\J(<Pk+l){t) - J(<Pk)(t)\\ < A / ds||^ fc+ i(a) - ^ fc (s)|| , (4.11) 

J t 

the last passage depending on Ea. (|4.4j) . Eq.s (|4.11|) and (|4.9f) ^ imply (j4.9|) i. + i . o 



4.7 Lemma. For all k, k' and n G N, it is 

\\ W - <p k \\ < £ e A ^ ^j^^ , h := min(A:, k') ; (4.12) 
so, (ipk) is a Cauchy sequence. 
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Proof. To prove Eq. (|4.12jl . it suffices to consider the case k' > k (so that h = k). Writing 
fk' ~ l Pk = Ylj^k&j+l ~ Vi) and usin S Eq-dHHB we get 

ii „^ ^ A^h-tpy 

WW - «Pk\\ < s 2^ - } • ( 4 - 13 ) 

j=k 

On the other hand, for each £ > 0, it is Ylj=k^ j /j- ^ EjSfcf7i ! < ^ £ fe A ! ! with 
£ = A(ti — to) w e obtain Ea. (|4.12|) . implying \\ipi-i — (fk\\ — ► for (k, k') — > oo ( 2 ). o 

Proof of Prop J3.4L (c^fc) being a Cauchy sequence, lim^i—^oo := exists in C([to, h], F); 
</? belongs to 2) because this set is closed. By the continuity of J ', we have 

= lim JXwO = lim ¥>fc+i = ¥ ; ( 4 -14) 

KH^+OO ki— >+oo 

thus, </> solves the Volterra problem (|2.21j) . Finally, the inequality \\tp(t) — ip ap (t)\\ < R(t) 
for all t £ [to,ti] is ensured by the definition of D. o 



5 An elementary application of Prop J3T4l 

The results we are presenting in the forthcoming Prop.s 15. H 15.21 are essentially known 
(see, e.g., for the case U(t) = 1 and dimF finite), but their derivation as a subcase 
of Prop l3,4l is instructive: the main idea is to use the zero function as an approximate 
solution. 

5.1 Proposition. Consider the Volterra problem \2.21\) . where: 

i) U is a strongly continuous linear semigroup, with an estimator u(t) := Ue~ Bt (U > 1, 
B > 0). 

ii) V is continuous and Lipschitz at fixed time on the strict subsets of its open domain. It 
is B(0,p) x [to,T) C DomV for some p £ (0, +oo] ; T £ (to, +oo], and 

V(0,t)=0 forte[t Q ,T). (5.1) 

There is a continuous function £ : [0,p) x [to,T) — > [0, +oo), non decreasing in the first 
variable, such that 

\\V(f,t)\\ < £(\\f\\,t) for (f,t) e B(0,p) x [t ,T). (5.2) 
Hi) The control problem 

R(t) = U£(R(t),t) - BR(t) , R(to) = U\\f \\, (5.3) 
has a solution R E C ([to, ij\r)> [0, p)), for some t^ S (to, T). 

Then, the Volterra problem \2.21\) has a solution <p £ C([to, ijv)> F) and, for all t in this 
interval, 

11^)11 < R(t) ■ (5-4) 

2 Incidentally we note that l|4.12|l could be improved, but this is unnecessary: this estimate is needed 
only to infer the Cauchy property of the sequence. 
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Proof. We apply Prop.s I3~H YAM with <p ap (i) '■= for t E [fQ,fjv); the function t in item 
ii) is a growth estimator for V from the approximate solution. The datum and differential 
errors are d((p ap ) = —fo, e((p ap )(t) = 0, so they admit the estimators 5 := ||/o||, e(i) := 0. 
With these estimators, problem (|3.17|) takes the form ()5.3|) . o 

The symbol tjy adopted here for the right extreme of the domain of R is chosen for future 
convenience; it emphasizes the dependence of this object on the norm of the initial datum. 
In the time independent case £(r,t) = £(r), Ea. (|5.3|) can be solved by the quadrature 
formula 

/ TTP( \ R =t ~ tQ; (5 ' 5) 
Ju\\M\ U£(r) - Br 

let us write the explicit solution in a simple case. 



5.2 Proposition. Let the previous assumptions be satisfied with to = 0, 5 = +oo, 
T = +oo and £(r,t) = Pr p [P > 0, p > 1). Then, the problem h5.'J\) has the solution 
R E C 1 ([0, t N ), [0, +oo)) defined hereafter. It is 

{ +oo if PC/l/oir 1 < B, 

(pc/f n/oir 1 ) */ pup n/oir 1 > p, (5 - 6) 
= r-(i/p)i og( i-p/«) vo<^<«, 

v 7 [ l/u if B = < u; ' 

R(t) := J- for all t E [0, t N ), (5.8) 

[ 1 - (PC/fll/oIlP- 1 - B)E B ((p - l)t) }—i 

[ e Bu -l 

Eb{u):={ — B~ ?/jB>0 ' (5.9) 
[ u ifB = 0. 

TTie function R has the following features. If PC/ p ||/o|| p l < B, R is decreasing and 
R(t) -^Ofort^ +oo. IfPUPWMr 1 = B, R(t) = const. = U\\f \\. If PU p \\f \\ p - 1 > B, 
R is increasing and R(t) — > +oo for t — > tjj. 

Proof. Everything follows in an elementary way from ([5.5)1 . o 

Remarks, i) A map V as in the Example of page E] has the properties required by the 
previous Proposition, if the function t i— > P(t) appearing in Ea. H2.15|) is bounded on the 
interval [to,T) under consideration. In this case, the growth of V from zero admits the 
estimator £(r,t) := Pr p , with P the sup of the function 1 1— > P(t). 

ii) Obviously enough: if Pt/ p ||/o|| p l < B, the Volterra problem (|2.21j) has a solution ip 
defined for all t E [0, +oo), and \\<p(t)\\ < R(t) -c for t -> +oo. If P?7 p ||/oir _1 = P, we 
have again a solution defined on [0, +oo), and < P||/o|| for all t. If Pf/ p ||/o|| p_1 > 

P, we can grant existence of a solution (p at least until the time t N in Eq. (|5.fi|) . and the 
bound ||y(t)|| < R(t) with R diverging at t N ; the result for this case can be applied to 
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blow up problems, to get a lower bound on the time of explosion of the solution and an 

upper bound on its growth. 

Example. We consider the Banach space 

F := C Q (R) := {/ : R -» R | / is continuous, f(x) -> for x -> 00} ; (5.10) 

11/11 := sup |/(x)| for / G F. (5.11) 
xgR 

We define a linear semigroup U : t G [0, +00) — ► G £(F) setting 

(U(t)f)(x):=f(x + t) (5.12) 

for all / G F; £Y is strongly continuous, and ||ZY(t)/|| = ||/|| ( 3 ). The generator of U is the 
operator 

A:=^ iCqHR) CF^F , /^/ a (5.13) 

where Qj (R) is the space of the C 1 functions / : R — > R such that /, vanish at infinity. 
We also introduce the function 

V:F^F, V{f):=f p (p > 1 integer), (5.14) 

which can be seen as a t-independent case of the Example on page|IJ with T(fi, f p ) := 

of course [|0>(/i, ...,/ p )[| < 
We consider the Volterra problem (|2.21|) with to := 0, and an arbitrary initial datum 
/o G F; in the special case /o G Cq(R), the corresponding Cauchy problem is 

tp(t) = <p(t) x + if(tf , ^(0) = f , (5.15) 

involving a first order wave equation with polynomial nonlinear ity The results of the last 
Prop 15 .21 can be applied in this framework with U = 1, B = and P = 1. For any /q G F, 
this Proposition ensures existence of the solution tp from time to 



1 



(intending t N := +00 if / = 0), and gives the bound 



(5.16) 



< R(t) , R(t) := ^ (5.17) 



l-ip- 1)11/011^]^ 



for all t G [0,tjv)- In this case, the accuracy of the estimates in Prop 15.21 can be checked 
in a very direct way, because the maximal solution of (|2.21|) is known; this is given by 

f(t)(x) = fo(x + t) _ for t G [0,^) 3 (5.18) 



l-(p-l)fo(x + t)P-H]r>- 



1 



3 In fact, U can be extended to a linear group, also defined for t < 0, but we do not emphasize this 
aspect: our general framework is designed for time evolution in the future. 
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tf:=sup{i>0 (p- 1) sup fo(x) p ~ 1 t < lj. (5.19) 

xeR 

If /o G Cq(R), the above tp also fulfils the Cauchy problem (j5.15j) : this implies the full 
equivalence of the Volterra and Cauchy problems for such an /q, in spite of the fact that 
item ii) of Prop |2.7l does not apply to the non reflexive Banach space C*o(R). The following 
facts occur: 

i) for p odd, or p even and sup^ fo(x) = sup x |/o(x)|, ■& equals the time t N in Ea. ()5.16j) : 
thus, Prop 15.21 gives the best possible lower bound on the existence time of the maximal 
solution. 

ii) For p even and < sup^ fo(x) < sup x |/o(ac)|, it is +oo > $ > tjy. 
hi) For p even and sup x fo(x) < < sup x |/o(x)|, it is +oo = $ > tjy. 

The accuracy of the growth estimate (j5.17|) is easily analysed by comparison with Eq. (|5.18|) . 
We think that better results would arise in cases ii) hi) by suitably generalizing the theory 
of approximate solutions to the framework of ordered Banach spaces ^Sl ; this will be done 
elsewhere. 



6 Approximate solutions on finite-dimensional submanifolds 
of F. 

Let us discuss a general scheme to construct accurate approximate solutions, and apply 
to it Prop IM. 41 to get information on the exact solution; a typical realization of this scheme 
is the Galerkin method, discussed in the sequel. 

The framework. From now on: IA is a strongly continuous linear semigroup with gen- 
erator A and an estimator u(t) := Ue~ Bt (U > 1, B G R); V : DomV C F x R — ► F is 
continuous and Lipschitz on the strict subsets of its open domain; (/o,io) 6 DomV. 
Our idea is to construct an approximate solution ip av for the Volterra problem (|2.21|) . 
lying on a finite-dimensional (linear or nonlinear) submanifold of F; we assume the latter 
to be coordinatized by some real parameters a k , labelled by a finite set of indices /. More 
precisely, we consider an injective C 1 map 

Q : DomQ CR 7 ->F, a = {a k ) keI ^ Q{a) , (6.1) 

with open domain, such that the partial derivatives 

d k Q{a) = J^(a) G F (k £ I) (6.2) 

are linearly independent for all a G DomQ; we regard ImQ as a multidimensional surface 
in F. We also suppose that 

ImQ C DomA , ImQ x [t ,T) C DomV (6.3) 

for some T G (to, +oo] , and ask Q to be continuous as a map to DomA with the graph 
norm. The approximate solution we consider has the form 

<p ap (t) := Q{a(t)) , a( ) G C 1 ^, *i), DomQ) , t h-> a{t) (6.4) 
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where [to, t\) C [to, T), and o( ) is a function determined in the sequel. Clearly, the datum 
and differential errors of ip ap are 

d(cp ap ) := g(a(to)) - f , (6.5) 

e(<p„)(t) := d k g(a(t))a k (t) - AG(a{t)) - V(G(a(t)),t) (6.6) 

(here and in the sequel, we employ the familiar Einstein's summation convention on re- 
peated indices). We prescribe a( ) to fulfil a Cauchy problem 

a(t) = X(a(t),t) , a(to) = ao 

where ao is an initial datum, and X : DomQ C R/ — > R 7 a continuous vector field; the 

criteria to fix ao and X are discussed later. 

For convenience, for all a S DomQ, a S R 7 , t £ [to,T) we put 

8(a) := \\G(a)- f \\ , e(a, a, t) := \\AG(a) + 7>(0(a), t) - d k g(a)d k || (6.7) 

i(a,t) := i(a,X(a,t),t) ; (6.8) 
then, the approximate solution <p ap admits the datum and differential error estimators 

5:=6(a ), e(t) := e(a(t),t) . (6.9) 

To conclude, we assume there are p £ (0, +oo] and a continuous function 

1 : [0, p) x Domg x [t , T) -> [0, +oo) , (r, a, t) ■-> l(r, a, t) , (6.10) 

non decreasing in the variable r, such that a £ Domg, t S [£ojT) and ||/ — <7(a)|| < p 
imply (/, i) S DomP and 

<i(||/-a(«)||,a,i) • (6.11) 

Then, the function 

t(r,t) :=£(r,a(t),t) (6.12) 
is a growth estimator for V from </? aJ) . The application of Prop.s l3~H 13.61 gives 

6.1 Proposition. Consider the equations 

a(t) = X{a(t),t) , a(t )=a , (6.13) 

= Ue(a(t),t) + Ul(R(t),a(t),t) - BR(t) , R(t ) = U8(a ) , (6.14) 

for the unknowns a( ) £ C ([io, ii), Domg), R G C 1 ([to, [0, p)). // (a( ),i?) zs a 
solution on some interval [to,ti) a^ic? ^ av (t) '■= g(a(t)), then the Volterra problem \2. 21]) 
has a solution ip on \t,Q,t\), and \\<p(t) — <p ap (t)\\ < R(t) on the same interval. o 
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Let us pass to the criteria for choosing X and ao- One of the most familiar is the Galerkin 
criterion (see, e.g., [I] or |15j): we will concentrate on it and will not discuss other ap- 
proaches (such as the variational methods often used for the Lagrangian or Hamiltonian 
evolution equations, see, e.g., [H]). The Galerkin choice for ao and X is the one minimizing 
the norms of the datum error and of the differential error (at any time): 

6.2 Definition. The vector field X and the datum ao fulfil the Galerkin criterion if 

e(a, d, t) = min! for a = X(a,t), (6.15) 

5(a) = min! for a = ao (6.16) 
(the symbol min! indicating the absolute minimum). o 

(Of course, condition (|6.16|) is trivially satisfied if /o = Giflo); then the absolute minimum 
of 5, attained at this point, is zero). 

Both equations (|6.15|) H6.16|) can be studied in a systematic way if F is a Hilbert space, 
say real, with an inner product < | > yielding the norm ||/|| := \J < f\f >; from now on 
we stick to this assumption. 
For all a G DomQ we introduce the matrix 

Qkl (a) :=<d k g(a)\d l g(a) > (k, lei), (6.17) 

which is symmetric and positive defined (recall the linear independence of the vectors 
dkQ(a)). As customary in tensor calculus, we denote the inverse matrix with g kh (a) 
and introduce the convention of "raising and lowering indices" with these matrices. In 
connection with this, it is worthy to write down the identities 

<d k g(a)\d h g(a) >=s\ ; (e.is) 

v k (a)d k g(a) = Vk {a)d k g{a) , 

< v k {a)d k g{a) I w h {a)d h g(a) >= v k (a)w k (a) = v k {a)w k {a) . (6.19) 

Here: d k G(a) := Q kh (a)dhG(a) G F; v k (a) (k G I) is a family of real numbers, v k (a) := 
Qkh{o)v h (a) and w k (a), w k (a) have a similar meaning. We apply all these notations to 
the discussion of the minimum problems (|6.15j) , (|6.16j) ; the solutions are given by the two 
forthcoming Propositions. 

6.3 Proposition. For any fixed (a,t) G Domg x [to,T), the function a i— > i 2 (a,d,t) is 
quadratic; it has a unique point of absolute minimum at 

d k = X k (a,t), X k (a,t):=<d k g(a)\Ag(a)+V(g(a),t)> . (6.20) 

The absolute minimum e(a, X(a,t),t) := e(a,t) is given by 

e(a,t) 2 = \\Ag(a)\\ 2 - < AG (a) | d k g(a) >< d k g(a) | AG (a) > + 

+ 2 < AG {a) | V(G(a),t) > -2 < AG (a) | d k G(a) >< d k G{a) | V(G(a),i) > + (6.21) 
+ \\V{G{a),t)\\ 2 - <V{G{a),t)\d k G{a) >< d k G(a) \V{G{a),t) > . 
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Proof. The general theory of Hilbert spaces tells us that, given a vector g and a closed 
vector subspace T of F, the problem 

find / in T such that \\g — f\\ = min! (6.22) 

has the unique solution 

/ = Tig; \\g — f\\ 2 = \\g\\ 2 — < g \ Tig >, II := the orthogonal projection of F — > T. (6.23) 

For any fixed (a,t), the problem we have in mind has just the form 1)6.22(1 : in this case 

T = T(o) := {h k d k g(a)) \ a G H 1 } C F , (6.24) 

(which represents the tangent subspace at G{a) of ImG), the unknown / is written 
X k (a,t)d k G(a) and 

g = AG(a)+V(G(a),t), U = U(a) =< d k G(a) \ ■ > d k G(a) ; (6.25) 

computing from (|6.23|) the solution / and the quantity \\g — f\\ 2 = e(a, t) 2 , we obtain Eq.s 
(IfHnl) (E23J. o 

6.4 Proposition. Assume that: 

i) Q is C 2 , DomG is convex; 

ii) the matrix Qm(cl)+ < Q{a) — /o | d^Gia) > is semipositive for all a G DomG (where d\ x 
are the second partial derivatives w.r.t. a k ,a l ); 

Hi) there is a point ao such that 

<d k G(a )\G(a )>=<d k G(a )\fo> . (6.26) 
Then, the absolute minimum of 5(a) is attained at a = a^. 

Proof. Everything follows computing the first and second derivatives of the function 
a — ► 8 2 (a) in Ea, (|6.7|) . Ea. (|6.26|) is the vanishing condition for the first derivatives, the 
other assumptions ensure that the stationary point a$ is of absolute minimum. o 

Remark. The geometrical meaning of Ea. (|6.26[) is Ti(ao)G(ao) = IT(ao)/o, where II(ao) 
is the orthogonal projection onto T(ao), see (|6.25|) . o 

The classical Galerkin method. All the previous formulas become very simple in the 
"classical" realization, where 

g : R 7 ^F, a^G(a)=a k e k , (6.27) 

(e k )kel linearly independent vectors of DomA, (e k ,t) G DomV for k G /, t G [to,T) . 
In this case ImG is a vector subspace, and we have the identities 

d k G(a) = const. = e k , g k i(a) = const. := Q ki , d k G(a) = const. = g kh e h := e k . (6.28) 

(Also, the tangent space at any point G(a) is T(a) = const. = ImG] further simplifications 
occur in the orthonormal case where Q k i = 5 k i, e k = e k ). 
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6.5 Proposition. The absolute minimum point for 5 and its value at this point are given 
by 

a k =<e k \f > , S(a ) = ||/ -age fc || . (6.29) 
(In particular, 5(ao) = if fo is in the linear subspace spanned by the family (e k )). 

Proof. Elementary (recall Eq. (|6.26j) ). o 

6.6 Proposition. Suppose V to be determined by a p-linear map 7 1 , as in the Example 
of page^Q Then, Eq.s \6.2U\) and \6.21\) become 

X k (a,t) =< e k \Aei > a l + < e k \ 9(e h , e lp , t) > a h ...a lp ; (6.30) 
e(a, t) 2 = ( < Aej \ Aei > - < Aej \ e k >< e k \ Aei > ^Ja j a l + 

+ 2( < Ae j \T(ei 1 ,...,ei p ,t) > - < Aej\e k >< e k \'S>{e h e ip , t) > ^a j a h ...a lp + (6.31) 
+ ( < ^(e jl ,-,e jp ,t)\?(e h ,...,ei p ,t) > + 

- < y(e h ,...,e jv ,t)\e k >< e k | 7{e h , e ip , t) > )a* ...a?*>a h ...a 1 " . 

In the r.h.s. of the last equation, the coefficients of a^a 1 and a^a ll ...a lp are zero if the 
subspace spanned by the family (e k ) is invariant under A (which occurs, in particular, if 
each e k is an eigenvector of A). 

Proof. Both Eq.s (j6.3QI6.31j) follow easily from d k Q{a) = e k and from the multilinearity 
of (P. In the second equation the coefficients of a^a 1 and a^a ll ...a lp are, respectively, 

< Aej | Aei > - < Aej | e k >< e k \ Aei >=< Aej \ Aei > - < Aej \ TIAei > , (6.32) 

< Aej\y(e h ,...,ei p ,t) > - < Mj\e k >< e k \ 9(e tl , ...,e lp ,t) >= 

=< Aej\T(e h ,...,ei p ,t) > - < UAej \ 7{e h , ...,ei p ,t) > , 

where II is the projection on the linear subspace spanned by (e k ). If this subspace is left 
invariant by A we have II^4e; = Aei for each I £ I, so the above coefficients are zero. o 

6.7 Proposition. Assume again V to be as in the Example of page^ Then, the growth 
of the function V starting from a point G(a) = a k e k admits this estimate, for all f 6 F 
and t G A; 

\\V(f,t) -V(g(a),t)\\ < i(\\f - g(a)\\,a,t) (6.33) 

l: [0,+oo) x R 1 x A -» [0,+oo), l(r, a, t) := P{t) ^ ( P \ y/a^aF P ~V . (6.34) 

j=i \ / 
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Proof. Apply Ea. (j2~T7|l with /' = Q(a) and i' = t; note that ||/'|| = y/a k a k . o 

In applications of the classical Galerkin method, some choices for the family (e^) have a 
special consideration. Apart from systems of eigenvectors of A, other choices occur in finite 
elements methods, which are strictly related to the idea of approximating the evolutionary 
problem by space discretization; in this case, (e k ) is typically a family of piecewise linear 
(or polynomial) "chapeau functions" related to some spatial grid, see e.g. [I]. 



7 Applications to the nonlinear heat equation. 

Our aim is to discuss the nonlinear heat equation / = f xx + f p with Dirichlet boundary 
conditions, for x ranging in (0, 7r) (we work in one space dimension only for simplicity). 
Let us introduce this equation in the framework of Sobolev spaces. 

Notations for Sobolev spaces. All functions on (0, tt) are real-valued; J-(0, vr) means 
^((0, 7r), R) for each functional class T. We consider the Hilbert space L 2 (0,vr), with the 
standard inner product < / 1 g >^2:= J™ dxf(x)g(x); here the functions 

s k (x) := ^ sm(kx) (fc € {1,2,3, ...» (7.1) 

form a complete orthonormal system. We introduce the Sobolev space 

^(0, vr) := {/ G L 2 (0, tt) I f x G L 2 (0, vr)} C C([0, vr]) , (7.2) 
f x denoting the distributional derivative of /; this is a Hilbert space with the inner product 

< / I 9 >■=< f\g>L* + < fx\ 9x >L2 , (7.3) 

yielding the norm ||/|| := y/ < f\f>. The inclusion indicated in (|7.2|) is well known, and 
allows to define /(0), f(ir) for all / G -ff 1 (0, tt); we fix the attention on the closed subspace 

F := Ht(0,ir) := {/ G H\0,tt) \ /(0) = /(vr) = 0} , (7.4) 

and equip it with the restriction of the inner product l|7.3|l . It turns out that 

oo 

F = {/GL 2 (0,vr) | ^A: 2 < Sfc |/>| 2 <+oo} . (7.5) 
fc=i 

The functions s k form a complete orthogonal system for this space: it is 

< s k | si >= (1 + k 2 )5 kU < s k \f >= (1 + k 2 ) < s k \f > L 2 V/ G F , 

oo 

< f\g>=Y,( 1 + k2 ) < f\ s k>L^<s k \g>^ V/, 9 gF. (7.6) 

k=l 



21 



Both H 1 (0, 7r) and F are known to be Banach algebras with respect to the pointwise 
product. We will use systematically the inequality (almost optimal, see the Appendix 

\\fg\\ < H/ll \\g\\ V/../6F. (7.7) 
The operator A.. This is the linear map 

A:=-fz on DomA := {/ G F | f xx G F} . (7.8) 

Of course, for all k, 

As k = -k 2 s k . (7.9) 
The operator A generates the strongly continuous linear semigroup IA on F, defined by 

oo 

U{t)f :=5>~ fc2 ' <s k \f> L 2 Sk for t G [0, +oo), / G F; (7.10) 
fc=i 

the above series is in fact convergent in F, and 



\Mt)f\\ 



E e- 2fc2t (l + fe2) < Sfc | / >2 2 < e -* ll/H . (7.11) 



fc=i 



The nonlinear function "P. This is defined by 

V : F -» F , := / p (p > 1 integer). (7.12) 

It belongs to the class of maps in the Example of page |1J and corresponds to the time 
independent p-linear map 

?:xPF^F, 7{f u ...J p ):=h...f p . (7.13) 

Of course, Eq.jZU) implies ^(/i, .., / p )|| < ||/i||...||/ p ||. 

The Volterra and Cauchy problems. We consider the Volterra problem (|2.21j) on 
F, with U as before, V(f,t) := the above defined V(f), to := and some initial datum 
/o G F; this reads 

<p(t) = U(t)f + I dsUit- s)p(s) p . (7.14) 
Jo 

From now on, we will denote with 

^M^F (7.15) 
the maximal solution. If fo G DomA, (|2.21l) is fully equivalent to the Cauchy problem 

<p{t) = <p{t) xx + ip(t) p , ^(0) = f (7.16) 
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(because F is reflexive and V Lipschitz on the strict subsets of its domain). 

Kaplan's blow up criterion. We specialize to the present framework a general citerion 
of Kaplan [5] for the blow up in a finite time of the solution of a nonlinear parabolic 
equation. To this purpose, we introduce the function 

Q:L 2 (0,vr)^R / h+ Q(/) := ± < sin \f > L 2= \ f dx sinx f(x) . (7.17) 

1 1 Jo 



7.1 Proposition. Consider the Volterra problem {7.1$ ; if 

fo > , Q(fo) > 1 , (7.18) 

then 

Proof. It is sketched in the Appendix^ adapting Kaplan's general argument. o 

From now on, t K will be referred to as the Kaplan time for the datum fo. 

Basic estimates on the solution. As a first step in our analysis, let us apply Prop.s 
15.11 and 15.21 to the Volterra problem. In the case we are considering, the semigroup has an 
estimator u(t) = Ue~ Bt with U = 1, B = 1; also, it is ||7>(/)|| < £(\\f\\) with £(r) := r p . 
Therefore, we have 

7.2 Proposition. For any initial datum fo ef, it is d > t N and \\ip(t)\\ < R(t) for all 
t G [0,t N ), with t N and R depending on the norm of fo in the following way: 

if ll/oll < 1, 

lo g (l-^) •/ Il/Oll >1, (7 ' 20) 

ll/oll 




[l-(||/ ||p-i-l)(e(p-i)«-l)]5=T 



If ll/oll < 1, R is decreasing and R(t) — > for t — > +oo. If \\fo\\ = 1, R(t) = 1 for all t. If 
|| /o|| > 1, R is increasing and R(t) — > +oo for t —* t^. o 

Summary of the previous results on An example. We have 

^ < for all / € F ; d < t K if / > 0, Q(/„) > 1, (7.22) 
with t N as in (|7.2U|) . t K as in ()7.19|) . Let us consider, in particular, the initial datum 

fo{x) := As%(x) = J^Asinx , (A > 0) ; (7.23) 
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then /o £ DomA, so we have a full equivalence of Q7.14|) with the Cauchy problem (|7,16|1 , 
It turns out that 

ll/oll C N :=^ = 0.7071.. ; Q(f ) = A C K := 2 J? = 1.595.. . (7.24) 

Cat 2 Ojf V vr 

Therefore, Eg. 1)7.22(1 with this choice of the datum tells us that 
6 = +oo if < A < C N ; t N <tf if C N < A < C K ; t N <$<t K if A > C K , (7.25) 

*« ■= lo S ( X " a^i) -A^+oc ^ §=r , (7.26) 

** := "iPT log ( X " l£) — ^i % ' (7 ' 27) 

It should be noted that (|7.22|) does not allow to establish whether i? is finite or infinite, for 
A in the interval (C N , C K }. In the rest of the Section, we will infer more precise estimates 
about $ by means of the Galerkin method, and also rediscuss its behaviour for large A. 

A Galerkin approach to the nonlinear heat equation. We apply the scheme of 
Sect El with 

I a finite subset of {1, 2, 3, ....} , Q{a) := a k s k for all a = (a k ) £ R/ (7.28) 

and s k the functions (|7.1|) . We refer, in particular, to the description given in the previous 
Section for the "classical" Galerkin method, to be employed with e k := s k , < \ > the inner 
product (|T3|) on F := H^(0,tt), and 

Qkl = (1 + k 2 ) S kl , g kl = 5 kl ; (7.29) 

these matrices are used to raise and lower indices. The vector field X, the error function 
e and the growth estimator function £ of Eq.s (|6.30|) (|6.31|) (|6.34j) are time independent, 
and given by 

X k (a) := -k 2 a k + < s k \s h ...si p > a h ...a l ? , (7.30) 
e(a) 2 := (< Sjl ...s jp \s h ...s h > - < Sjl ...s jp \8 h >< s k \s h ...s lp >)a jl ...a jp a h ...a lp , (7.31) 

£(r, a) := £ ( P ) P ~V (7.32) 

(see the observation following Ea. (|7.13|) ; also, note that a k a k = ^2 k€l (l + k 2 ){a k ) 2 ). This 
amount of information (completed by Eq. (|6.29|) for the initial datum) must be inserted into 
the general scheme of Prop 16. ll the solution of the finite dimensional system (|6.13|) . (|6.14|) 
appearing therein provides simultaneously: 

i) a pair of functions a( ), R{ ) on an interval [0, t G ), the former giving the approximate 
solution <p a p(t) '■= a k (t)s k . In the sequel t G will be called the Galerkin time; 
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ii) an assurance that the Volterra problem l|7.14|) has a maximal solution ip on [0, 9) D 
[0, t G ), and a bound \\<p(t) — (p ap (t)\\ < R(t) for all t < t G . 

Introducing an example. We assume 

p := 2 , f as in QFZty . (7.33) 

Problem (1 7, 14117,16]) will be treated with a "two-modes" application of the Galerkin method; 
more precisely, we will work on the linear submanifold spanned by (sfc)fcgj, setting 

I:={1,3}, a:=a\ 7 := a 3 ; (7.34) 

(the choice J = {1, 2, 3} would not yield any improvement, because the function 1 1—* a 2 (t) 
would be ultimately found to be zero). The vector field, the error function and the growth 
function of Eq.s (17301 (T73H (|732l are given by 

^(a,7) = -a + /5(^ 2 -^7 + |7 2 ) , (7.35) 
2"/ 8 „ 144 



X 7 (a,7) := -97 + \/ I — -a 2 H «7 + ~7 2 

w 7T 3 V 15 35 ' 9 ' 



e>,7) 2 = hr- " « 4 + TT^T " — « 3 7+ (7-36) 



i2 ^ 7 512 \ 4 / 34816 10 
2k ~ 15vr 3 / a + V315vr 3 ~ "tt 

/46 12172288 \ 2 2 22528 3 /39 3247616 \ 4 
+ V~ ~ 33075vr 3 J ° 7 175vr 3 ° 7 + V2^ ~ 99225vr 3 J 7 ' 

7, r) = r 2 + 2\/2a 2 + IO7 2 r . (7.37) 
According to (|6.29j) . the initial conditions for a(t) and j(t) are, respectively, 

< s 1 I /o >= A , < s 3 I /o >= ; (7.38) 

the corresponding datum error is zero. In conclusion, we have to study the system 

a = X a (a,j), 7 = X 7 (a, 7 ), R = e(a, 7) + £(a, 7, R) - R , (7.39) 

a(0) = A , 7(0) = , R(0) = , 

for the unknown functions t 1— > a(t),j(t),R(t). This cannot be solved analytically, but 
can be easily treated by any package for the numerical solution of ordinary differential 
equations; an integration algorithm with adaptative control of the step size gives an excel- 
lent approximation for the solution of (|7,39|) , also including the evaluation of its existence 
time. All statements that follow about the system ()7,39|) are based on the MATHEMAT- 
ICA package; thus, expression such as "the solution of (|7.39j) " ■ etc., always indicate the 
MATHEMATICA output (of which we report the first digits). 

New estimates on the existence time i9. We have the bounds 

t a < ■& for all A > 0, ■& < t K for A > C K (7.40) 
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involving the Galerkin and Kaplan times, defined previously (the latter depends on C h 
1.595..). It is found that 



t G = +00 (whence 1? = +00) for < A < Q G , Q G = 1.056.. ; 
a(t),j(t),R(t) for t -> +00 and < A < C G . 



(7.41) 
(7.42) 



It should be noted that the bound C G on A for the global existence of (|7.16|) is better than 
the previously derived bound C N = 0.7071... For larger values of A, the existence time 
t G for (|7.39[) is finite. The forthcoming Table reports t G for some values of A above the 
Kaplan critical value C K , as well as the corresponding values of t K . It also reports 



77 := 



t K + t c 



(7.43) 



which is the relative uncertainty about the actual existence time 1? of (|7,16|) corresponding 
to the lower and upper bounds t G , t K . 



A 


to 


t K 


V 


< 1.056.. 


+00 






1.60 


1.104.. 


5.935.. 


0.6861.. 


2 


0.7730.. 


1.598.. 


0.3481.. 


4 


0.3138.. 


0.5090.. 


0.2372.. 


10 


0.1112.. 


0.1738.. 


0.2196.. 


20 


0.05340.. 


0.08315.. 


0.2177.. 



(7.44) 



The A — > +00 limit for the previous estimates is easily discussed. To determine the 
asymptotics of t G , we reexpress the unknown functions a(t), 7(t) and R(t) in terms of 
three rescaled functions t — ► a(t), c(t),R(t), depending on t := At and defined by 



a(t) = A a(At) , -y(t) = A c(At) , 
Then, the system ()7.39() becomes (with ' := d/dt) 



R(t) = A R(At) . 



(7.45) 



A 



+ 



16 



72 



— ac H c" 

15 35 



9c 
~A 



+ 



15 



:a 2 + 



144 8 , 

ac H — c x 

35 9 



R' = e(a, c) +£(a, c,R) 
a(0) = 1 , c(0) = , 



R 

~ A 

R(0) = . 



(7.46) 



In the A — ► +00 limit, the terms a/A, c/A and R/A can be neglected and the outcoming 
system ()7.46|L C is A- independent. The numerical treatment of this limit system shows 
that the solution t — > a(t), c(t), R(t) exists for t G [0, C G ), where 



1.026.. 



(7.47) 
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Returning to the standard time t = t/A, we conclude that 



to ~ -f- for A -» +00 . (7.48) 

It should be noted that all values of ^4 in the previous Table are seen empirically to fulfil 
the inequality 

with i G very close to the r.h.s. To compare these results with the Kaplan upper bound 
recall that, for all A > C K = 1.596.., 

t K = - log (l - tf\ ~a-, + oo • (7.50) 

Due to Q7.48|) (|7.50[> . the relative uncertainty Q7.43|) has the limit 

C — C 

rj -» — = 0.2173.. for A +00. (7.51) 
+ C G 

As a matter of fact, in this limit case one can find directly the asymptotics for the actual 
existence time i? of (|7.14I7.16|) . In fact, if one writes the maximal solution ip as 

ip(t) = A X (At) (7.52) 

one obtains for x the Cauchy problem 

X '(t) = j X (t)xx + X 2 (t) , x (0)(x) = y|sinx . (7.53) 

For A — ► +00, the differential equation becomes x'fa) = X 2 (t)> and the solution is 

X(t)(x) = - ^^ n (") for t e [0, y^) . ( 7 .54) 

1 — y/Z/nt smx 

so, returning to the standard time t = t/A, we conclude 



1.253.. 

1? ~A^+oo — ^— = A ■ (7.55) 



The constant \/tt/2 is fairly close to the arithmetic mean of the costants C G and C K , 
which appear in the asymptotic expressions (|7.48j) . (|7.50|) for the lower and upper bounds 
t G , t K . Thus, the actual existence time •& should be close to the arithmetic mean of t G and 
t K if A is sufficiently large. An attack to the Cauchy problem 1)7. 16j) that we performed 
approximating d 2 /dx 2 by finite differences seems to indicate that this actually occurs for 
all A > 2. 
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Analysis of the Galerkin solution. After solving the system (|7.39j) for a given A, one 
constructs the corresponding approximate solution for (|7.16|) such that, for t E [0,i G ), 





<fa P {t)(x) =a(t)A/-sinx + 7(t)W-sin(3x); ||<M*)II = VM*) 2 + Wjffi ■ (7-56) 



The system (|7.39j) also gives a function R such that ||(/? ap (i) — (p(t)\\ < R(t) for t in the 
same interval; let us illustrate the behaviour of the above functions for two values of A. 
Case A — 1. It is t G = +00, which implies 1? = +00. Figures ^ 121 HI give the graphs of 
the functions 1 1— > a(t), ~y(t), \\(p ap (t)\\ and R(t), all converging to zero for t — ► +00. Figure 
131 gives the function x G (0,7r) i— > <^ ap (t)(x) at three fixed times. For the exact solution ip 
of (|7.16() . we infer 



FigureEJis a graph of the relative bound i?(t)/||<^ ap (t)|| in a time interval where it is fairly 
little. 

Case A = 4. The Galerkin system (|7.39J) has a finite existence time t G = 0.3138... Figures 
l^l llUI give information of the same kind as the figures of the case A = 1, but describe a 
qualitatively different behaviour; in particular, the function 1 1— * R(t) diverges for t — > t G . 



Mt)\\ <IIMt)ll+fl(*)-o 



for t 



+00 . 



(7.57) 
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1 2 3 4 5 

Figure 1: A = l. Graph of a(t). 




Figure 2: A = 1. Graph of j(t). 




Figure 3: A = 1. Graphs of ip ap (t)(x) for x € (0, 7r) and £ = 1 (continuous line), t = 2 
(short dashes), t = 3 (long dashes). 





Figure 4: ^4=1. Graphs of ||(/9 ap (t)|| (contin- 
uous line) and R(t) (dashed line). 



0.25 0.5 0.75 1 1.25 1.5 

Figure 5: A = l. Graph of R(t)/\\<p ap (t)\\. 
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Figure 6: A = 4. Graph of a(t). 
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Figure 7: A = 4. Graph of -y(t). 
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Figure 8: A = 4. Graphs of ip ap (t)(x) for x G (0,7r) and t = 0.1 (continuous line), t = 0.2 
(short dashes), t = 0.3 (long dashes). 



0.05 0.1 0.15 0.2 0.25 



Figure 9: A = 4. Graphs of ||<^ op (£)|| (contin- 
uous line) and R(t) (dashed line). 
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Figure 10: A = 4. Graph of R(t)/\\ip ap (t)\\. 
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A Appendix. Sobolev spaces and pointwise product. 

We consider the space -ff 1 (R, C) := {/ : R — » C | /, / x G L 2 (R)} with the inner product 
< / 1 9 >■=< f | 9 >l 2 + < fx I 9x >l 2 an d the corresponding norm || ||. This is useful to 
treat the space Hq(0, tt) of Sect [7] (made of reai functions): in fact, there is an R-linear, 
norm preserving inclusion 

flj(0,7r) <ZH\R,C) , (A.l) 

where each / G Hq(0,it) is extended to the full real axis setting f(x) := for x G" (0,7r). 

Both Hq(0,it) and i/ 1 (R, C) are closed under the pointwise product, and < const. 
Il/H ||<?|| for all functions therein [TI|. We claim the following: 

A.l Proposition. Consider the sharp (i.e., the minimum) constants L,M in the in- 
equalities 

II/5II < L\\f\\ \\g\\ for all f,g€ H^ir) ; (A.2) 

||/<?ll < M\\f\\ \\g\\ for all f,g£ H\R, C) . (A.3) 

Then 

0.811 < L < M < 1 . (A.4) 
Proof, i) L < M follows readily from (jA.l|) . 

ii) A lower bound for L follows applying (|A.2|) with f = g = fx, where 

fx(x) := e-AI^-TT/21 _ e -A^/2 (A > 0) , (A.5) 

The norms ||/a||, ||/?|| are computed in an elementary way, and we get a minorant of L 
for each A. The best lower bound is attained for A close to 1.55, and implies L > 0.811. 

iii) Let us prove that M < 1. To this purpose, we employ for the complex functions / on 
R the Fourier transform (J : f)(k) = Jj^dx e~ tkx f(x); this yields the representation 

H\R, C) = {/ G L 2 (R, C) | +oo > f dk{l + k 2 )\Ff(k)\ 2 = ||/|| 2 } , (A.6) 

and sends pointwise product into (1/\/2tt)x the convolution product * . Consider any two 
functions /, g G -ff x (R, C). Then the following holds: 

||/<7|| 2 =/ dk(l + k 2 )\f(fg)(k)\ 2 = ±- [ dk(l + k 2 )\(Ff*Fg)(k)\ 2 ; (A.7) 
(Ff * Fg){k) = [ dh Tf(k - h)Tg{h) = (A.8) 

Jn 

= [ dh — L_ ; (Jl + (k- h) 2 Tfik - h)y/l + h 2 Fg(h)) ; 

JR ^/i + (k-h) 2 VT+h 2 V v v ;v yv v 

±-\{Tf*Tg){k)\ 2 <C{k)P(k) , (A.9) 
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2vr 7r (1 + (A; - h) 2 ){\ + h 2 ) 4 + k 2 

P(k):=[ dh (l + (k-h) 2 )\Ff(k-h)\ 2 (l + h 2 )\Fg(h)\ 2 . (A.ll) 
./R 

Ea. (|A.9|) follows from Holder's inequality | J dh UV\ 2 < (J dh\U\ 2 ) (J dh \V\ 2 ); inserting 
1)A.9|) into Ea. (|A.7|) . we get 

||/#|| 2 < f dk(l + k 2 )C(k)P(k) < ( sup (1 + k 2 )C(k) ) f dk P(k) = (A.12) 
^R \*eR / -/R 

= ix||/|| 2 NI 2 - 

In we have discussed the constants for more general inequalities related to the pointwise 
product and to the spaces H n (R d ,C). The upper bound M < 1 derived now improves the 
result arising from jTJ in the special case of the inequality (jA.3|) : the method employed 
here to bind M develops in a fully quantitative way an idea suggested in |13| . 

B Appendix. Proof of Prop JTTTl 

We keep all notations of Sect|7| The proof consists of the following steps: 

i) The function Q of Ea. (|7.17j) can be seen as a continuous linear form, both on L 2 (0,7r) 
and on F. For all / G F and t G [0, +oo), we easily infer from (|7,1U|) that 

Q(U(t)f) = e-tQtf) . (B.l) 

ii) Let / G L 2 p{0,tt)(c L 2 (0,tt)) and / > 0; then 

Q(f) < Q(f p ) 1/P ■ (B.2) 

This follows taking q such that 1/p + 1/q = 1, and writing Q(f) = L dx u(x)v(x) with 
u(x) := (^sinx) 1 / 9 , v(x) := (| sinx) 1 / p /(x); Holder's inequality J uv < (J v p ) l l p 
yields Ea.(lB~2l) ( 4 ). 

iii) We consider the maximal solution 99 : [0, 9) — * F of the Volterra problem (|7.14l) with 
datum /o G F, assuming /o > and Q(fo) > 1; the nonnegativity of /o implies ^c(t) > 
for all £ (see, e.g., [2]). We define the (continuous) function 

t G [0, 9) \-* Q(t) := Q(<p(t)) , (B.3) 

and note that 

Q(t) > , Q(t) > e-*Q(/ ) + /* e -(*- s )Q(s) p ; (B.4) 
Jo 

4 Ea. llB.2fl is optimal, in this sense: the best constant in the inequality Q(f) < CQ(f p ) 1 ^ p for all 
nonnegative L 2p functions is C = 1. This is true even if we restrict the inequality to much smaller classes, 
such as the nonnegative C°° , compactly supported functions / on (0, tt) . 
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the first bound follows from (p(t) > 0, and the second one from (|7.14j) (|B,1|) (|B.2|) . 

iv) For each n G N, we define a continuous function S n : [0, i?) — > R by 

Sb(t) := e-*Q(f ) , S„+i(t) := e -*Q(/ ) + / ds e^^S^f ; (B.5) 

j o 

it is proved recursively that, for all n £ N and t G [0, i?), 

> S n (t) > (B.6) 

(the first inequality depends on (|B.4|) . the second one is elementary). It is easily checked 
that the sequence of functions (S n ) is a Cauchy sequence in the topology of the uniform 
convergence on all compact subintervals [0, r] C [0,19); its n — * +oo limit is a continuous 
function S : [0, i?) — > R such that, for all t in this interval, 

S{t) = e-*Q(/ ) + T (is e-(*-^(s) p , Q(t) > 5(t) > 0. (B.7) 

From the above integral equation, we see that S is in fact C , and fulfils the Cauchy 
problem 

S(t) = S(t) (S(ty- 1 - l) , 5(0) = Q(/ ) ; (B.8) 

this has a unique maximal solution (for nonnegative times), denoted again with 1 1— > S(t), 
which extends the function considered up to now on [0, and is given by 

S(t) j r r+oo ^ r 

t fortG[0,t K ), t K := I - 7 _ r — — ; (B.9) 



Wo) 1 " !) " " " ' " ' ^Q(/o) rir '' 

furthermore, 5(t) — ► +cxd for t —* tg. Computing the last integral, we see that t K has the 
expression 1)7. 19|) in the statement of the theorem; we know that i? < t K , so the proof of 
Prop 17. ll is concluded. o 
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